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Abstract 

We present a brief overview of integrability of nonlinear ordinary 
and partial differential equations with a focus on the Painleve prop- 
erty: an ODE of second order has the Painleve property if the only 
movable singularities connected to this equation are single poles. The 
importance of this property can be seen from the Ablowitz-Ramani- 
Segur conhecture that states that a nonlinear PDE is solvable by in- 
verse scattering transformation only if each nonlinear ODE obtained 
by exact reduction of this PDE has the Painleve property. The Painleve 
property motivated motivated much research on obtaining exact so- 
lutions on nonlinear PDEs and leaded in particular to the method of 
simplest equation. A version of this method called modified method 
of simplest equation is discussed below. 
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1 Several words on the integrability of differ- 
ential equations 

The development of the modern natural sciences leaded to wide use of non- 
linear models for natural and social phenomena PP-|1]. Thus the problems 
connected to obtaining exact solutions of the model nonlinear differential 
equations are very actual [5]- [9]. For many researchers the question about the 
integrability of nonlinear differential equations or systems of such equations 
is not very clear. They have probably heard is that there exist a complicated 
method called Inverse Scattering Transformation (1ST) (TUl [H] that allows 
obtaining exact solutions of some nonlinear PDEs. May be the researchers 
have heard that integration by 1ST can be performed if the corresponding 
nonlinear differential equation is somehow connected to a property called 
Painleve property that states something about the singularities of the solu- 
tion in the complex plane. Many researchers are puzzled by this fact as they 
have to solve differential equations which contain real variables. We shall 
discuss below the connection between the methodology of the test for the 
Painleve property and the research on direct methods for obtaining exact 
solutions of nonlinear differential equations. 

Given an differential equation or a system of equations how one can tell 
appriori whether or not they are integrable? For a second order ordinary 
differential equation the existence of a constant first integral (which for an 
example is the energy for the many problems connected to the physics) leads 
to a reduction of the order of the equation by 1. Thus the task for finding 
a solution is reduced to quadrature. For second order systems especially if 
they not conservative the integration is more difficult. Such systems possess 
solution but the problem for finding analytical representation of the solution 
is very difficult to solve. It is extremely useful to find integrals of motion 
connected to the solved equations. But how to find integrals of motion if 
they exist? The existence of first integrals may be connected to the analytical 
structure of the studied differential equation as it was implied by the research 
of Sofia Kovalevskaya. 

2 The research of Kovalevskaya and Painleve 

Kovalevskaya won the Bordin prize of the Paris Academy of Sciences in 1888 
for the solution of the system of six first-order nonlinear ODEs that describe 
the motion of a heavy top rotation about a fixed point. Kovalevskaya solved 
this mechanical problem by the methods of complex variable theory (she ob- 
tained conditions under which the only movable singularities of the solutions 
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in the complex plane are single poles). Movable singularity is a singularity 
whose location depends on the initial conditions. Let us consider the simple 
nonlinear equation in the complex plane [12] {z below denotes a complex 
number) 

(1) £+/'w-°- 

It has the solution 

(2) f{z) = ^; z, ^ 




z + z,' f{0) 

In other words f{z) has a singularity at ^ = 2:9 and this singularity is movable 
as it depends on the initial condition. Painleve continued this research on 
the connection between analytic structure and integrability by analysis of 
the class of second order differential equations 

(3) 

where F is analytic function in z and rational in y and Painleve obtained 
50 types of such equations whose only movable singularities were simple poles. 
44 of these equations can be integrated in terms of known functions. The 
remained six equations are called Painleve transcendents and they can not 
be integrated by quadratures. The simplest two Painleve transcendents are 

PI: = 6?/^ + x; P2: ^ = 2y^ + xy + tt; tt: free parameter . 

The solutions even for simple ODEs in complex plane are not known but 
the nature of their movable singularities can be determined by studying the 
local properties of the solutions. Let as above the function F is analytic in z 
and rational in other arguments. The behavior of the movable singularities 
of the n-th order ODE 

(4) ^^^/^ dy \ 
^ ' dz- [dz^-^'-'-'dz'^' ) 

can be determined by leading order analysis, i.e., by the ansatz 

(5) y{z) = a{z - zoT. 

Let us discuss an equation which contains as particular case the first Painleve 
transcendent {A below is a parameter) [12] 

(6) pi = 6y'' + Az. 



dz 
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We substitute Eq.([S]) in Eq.dH]) and obtain as a result a relationship that 
contains the powers a — 2; 2a;; and a. In order to obtain a non-zero expression 
of kind ([5]) we have to balance at least two of the powers. Taking into an 
account that the highest order singularity arises from the lowest (negative) 
value of a the appropriate balance is: a — 2 = 2a — )■ a = —2. Then at 
z close to zn the leading order of the solution is -, — i.e. second order 
(movable simple) pole. In order to characterize better the behavior in the 
neighborhood of the singularity a local expansion of the solution has to be 
considered. If the singularity is a movable pole the expansion is a Laurent 
series. In our case (the second-order pole above) the candidate for Laurent 
series is 

oo 

(7) y{z) = ^ai{z - zof'^ . 

1=0 

After the substitution of Eq.(I7]) in the equation ([6]) one obtains a system of 
recursion equations (ao = 1): caii + 1)(« — 6) = 6 X]l=i cti-i^i + Aaj_2- The 
solution of the system of equations is: ai = 0; 02 = c^s = 0; 04 = — 
05 = 0; O.ae = Vla^a/^^ Aa/^ = 0. Thus the series ([7]) have two free parameters 
Gq and Zq (the arbitrariness of the pole position). As the investigated ODE 
is of second order its general solution contains two free parameters. These free 
parameters are contained in the local series expansion, i.e., ([7]) are Laurent 
series of the general solution indeed. In addition in the neighborhood of the 
movable singularity zq the solution of Eq.([n]) behaves as second order pole, 
i.e., the equation has the Painleve property (which is a strong hint for its 
integrability) . 

The powers of (z — zq) at which the arbitrary coefficient appear are called 
resonances. At the resonances one obtains a compatibility condition which 
must be satisfied in order to ensure the arbitrariness of the corresponding 
coefficient. In our case the resonance occurred at z = 6 and the compatibility 
condition is: 120204 + Aa^ = 0. 

The leading order of the solution behavior near the pole as well as the 
value of the resonance can be integer number (and then the pole is simple 
and the equation has the Painleve property) but these two numbers can be 
also non- integer (for an example irrational or even complex). In this case the 
singularity is not a simple pole and the equation does not have the Painleve 
property (which is a hint of non-integrability from the point of view of the 
research of Kovalevskaya and Painleve) . 

In some cases the Painleve property can be established only for selected 
values of the parameters of the studied nonlinear ODEs. In such case one 
can attempt the \E'-series technique that will not be discussed here. 
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3 Weiss- Tabor-Carnevale research on Painleve 
property for nonlinear PDEs 

Weiss, Tabor, and Carnevale [13] continued the research of Painleve in the 
direction of finding hints for obtaining exact solutions of nonlinear partial dif- 
ferential equations. They work was motivated by the conjecture of Ablowitz, 
Ramani and Segur [13] that a nonlinear PDE is solvable by inverse scatter- 
ing transform (1ST) only if each nonlinear ODE obtained by exact reduction 
of this PDE has the Painleve property. In other words we have to test all 
reductions of the PDE to ODEs for Painleve property and if all ODEs have 
the Painleve property then the corresponding nonlinear PDE has to be in- 
tegrable by 1ST. But what are all these reductions? We don't know them. 
Then one has to follow another way and one possible way the the generalized 
Laurent expansion approach of Weiss, Tabor and Carnevale. The key point 
in this approach is that as difference to the case of complex function of one 
complex variable for the case of functions of several complex variables the 
singularities can't be isolated and these singularities are located in singular 
manifolds. If f{zi, . . . ,Zn) is a meromorphic function of the complex vari- 
ables zi, . . . ,Zn the singularities of / are located in analytic manifolds (called 
singular manifolds) of dimension 2n — 2 which are determined by conditions 
like 

(8) V(^l,---,^n) = 

where ip is analytic function in the neighborhood of the manifold defined 
by Eq.([H]). Weiss, Tabor and Carnevale have proposed generalized Laurent 
expansion of analytic function f{zi,..., Zn) as follows 

oo 

(9) f{zi,...,Zn) = —Y.Ui(t)' 

where ^(^i, . . . , z„) and Ui(zi, . . . , z„) are analytic functions of zi, . . . , 2„ in 
the neighborhood of the manifold given by Eq.®. a is an integer. Difference 
to the classic Laurent series is that instead oi z — Zq we have the function 
(j) and Ui are functions and not constant coefficients as in the classic case. 
Additional requirement to the function (p is that its gradients do not vanish 
on the singular manifold. 

The test for the generalized Painleve property happens as in the case 
of ordinary differential equations. First Eq.([9]) is substituted in the studied 
nonlinear PDE and the possible values of a and the recursion relations for the 
function Ui are determined. The recursion relations are coupled PDFs con- 
taining the functions and Ui. As in the case of ODEs discussed above there 
will be powers of for which the corresponding Ui are arbitrary (resonances). 
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If some Ui are not arbitrary at the resonance some kind of generalized \l/-series 
can be introduced to ensure their arbitrariness. 



4 Application to systems connected to fluid 
mechanics 

We shall discuss the application the methodology for detecting Painleve prop- 
erty to the Lorenz system and one of its generalizations, namely 

dx 

— = -ax + cry - uyz, 

10 — = rx — v — rnxz, — = —oz + xy. 

^ ' dt ^ ' dt 

where a, co, r, m, and b are parameters. When u = 0; m = 1 the system 
ffTOj) is reduced to the classical Lorenz system well known from meteorology, 
nonlinear dynamics and chaos theory. The classic Lorenz system was tested 
for presence of Painleve property in [15]. Let us transform the equations of 
the system flTU]) by setting x = ^] y = z = t = eT; e = The 
system (fTOj) becomes 

-aeX + Y- ^YZ. 



dT 



:il) — = X-eY - mXZ, — = -heZ + XY. 

' ' dT ' dT 



We shall investigate the system of equations flTTl) for presence of Painleve 
property. 

First of all we have to perform the leading order analysis by setting 
(12) ^ = j^'^ = T^' ^ " T7' 

where a, b, c and a, (3, 7 are parameters. The substitution of Eqs. f[T2|) in the 
system f lTTj) leads to the following relationships after balancing the largest 
powers in the denominators of the resulting three equations: a = /3 = 7 = 1. 
This result is different from the case of classical Lorenz system where m = 1 
and u = (the last relationships leads to significant changes of the right-hand 
side of the first equation from (fTTj) - the nonlinearity there is removed). For 
this case the substitution of Eqs. (|T2|) in the system (fTTj) leads to the result 
a = 1; /9 = 7 = 2. For the case of generalized Lorenz system we obtain 
the following system of equation for the parameters a, b, c (after equating 



6 



the numerators of the two largest powers from the denominators from the 
equations for the leading order analysis): 

(13) — a + = 0, — 6 + mac = 0, —c — ab = 0. 

The non-trivial solution of the system of equations f|T3l) is a = ±:^; & = 

y ' y mo; 

What follows are the power series expansions around the singularities. 
On the basis of the above results we shall consider the following expansions 
for the three unknown functions around the singularity of the solution which 
we assume to be at (Xq, Yq.Zq): 

V^(T - To) ^t^o ' v^(r - To) ^t^o °^ 

oo 



(14) 



These series are different from the series for the classical Lorenz system. For 
the classical Lorenz system the system of equations for the parameters a, 6, c 
is 

(15) a + 6 = 0, ac - 26 = 0, 2c + a6 = 0. 

which solution is a = ±2?; 6 = =p2z; c = —2, and the expansions for X, F, Z 
around the singularity become 

2? °° —2? °° 

U - -foj j=o - ^oj j=o 

r) oo 

The next step is obtaining of the recurrence equations. We have to sub- 
stitute the expansion f|T^ in the generalized Lorenz system ffTOj) and by 
equalizing to the coefficients of the powers (T — To)~", n = — 1, 0, 1, . . . we 
shall obtain a system of nonlinear recurrent equations for the coefficients aj, 
bj and Cj. Setting to the coefficients of (T — To)~^ in the resulting equations 
we obtain that ao = 6o = Co = 1. Setting to the coefficients of (T — To)~^ 
in the resulting equations we obtain the following system for the coefficients 
ai, hi, ci: 

a\fem{^ — a/cJ) -|- ci + 6i = 0, ot\C\\ a\a — y/ume = 0, 
(17) ai + bi + eb = 0. 
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The solution of this system is 



bi 

:i8) ci 



oj m 

— e + a\ a\ ae — eh 

V em V euj 



ae — a 



e + a 



1 

' 2 

+ (T 



a 
me 



\fujem — a 
— ae + eh 



me 



+ ae — he 



The situation with classic Lorenz system (a; = 0, m = 1) is different 
again. In this case we have to substitute the series fll6p in f llOp (with uj = 
and m = 1). Setting to the coefficient of (T — Tq)^^ in the first resulting 
equation and the coefficients of (T — To)~^ in the second ant third resulting 
equation we obtain the following result ao = ho = cq = 1. Next we have 
to set to the coefficient of (T — Tq)~^ in the first resulting equation and 
the coefficients of (T — Tq)~^ in the second ant third resulting equation. We 
obtain the following system of equations for ai, hi and Ci: 

2aei + 2ihi = 0, hi - e - 2ci - 2ai = 0, Ci - 2hi - 2ai - eh = 0. 

The solution of this system is 

(3cr-26-l)e , h - 1 - 3a 



(19) 



ai 



hi 



-ea, ci 



6 ' ^ ' ^ 3 

For j = 2, 3, . . . for the generalized Lorenz system we obtain the system of 
recurrence equations 



(j 



1 N 

\)aj H ^a,_i 



m 



ey/e 



(j - l)hj i=aj-i + —j=bj-i + 



ma' 

re 



m^/oj 



(20) 



1 / ■ . N be 1 
— (j - l)Cj + —Cj-l + 
m m \ m 



1=0 

j 

1=0 

j 

E (^ibj-i 

1=0 



0, 



0. 



Let us obtain the first relationship from the system (120|) in more detail. 
The substitution of the series (fT^ in the system of generalized Lorenz equa- 
tions ( 1TT|) leads to the relationship 



E(^-^o)^'' 



(21) 



j=0 

re 



m 



\l+k-2 



1=0 k=0 
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From Eq. fl?!]) we want to select the terms of the same power of (T — To). 
This means that we must have j — 2 = l + k — 2, i.e., j = l + k. In other words 
from the double sum in (l2Ti) we have to select the terms for which k = j — I 
and in addition we have to take into an account that A; > which leads to 
cutting the sum with respect to / from / = to / = j. Then the terms from 
the double sum that participate in the recurrence relationship are given by 
{6ij is the delta-symbol of Kronecker) 

j oo j 

(22) J2 E bicA,-iiT - ToY^"-' ^ E bic,-i{T - T,y-\ 

Z=0 fc=0 /=0 

We note that in the last relationship j was fixed. Putting all together for 
fixed value of j in Eq. fl2Tl) we obtain the first relationship from fl20l) . The 
remaining relationships from fl20|) are obtained in the same way. 

For comparison the situation in the classic Lorenz system is simpler. For 
j = 2, 3, ... we obtain the system of recurrence equations 

(j — l)aj + hj + aeaj^i = 0, 

i-i 

2aj + {j - 2)hj + 2cj + 2 ^ flj-fcCfc + aj-2 + thj-i = 0, 

k=l 

i-i 

(23) 2aj + 2bj + {j - 20cj + 2 ^ aj-kbk + ebcj-i = 0. 

fe=i 

For j = 2 and j = 4: there are resonances that lead to the following relation- 
ships among the parameters of the equations of the classic Lorenz system 

(6- l)[24 + 57(a 

(24) 6{1 + b-a)alc + e{2(j-b 

The last relationships fix conditions on the values of e, a and /3 for which one 
the classical Lorenz system is integrable. As e participates only in the third 
equations from fl24|) from the first two relationships there one determined a 
and b and then the last relationship leads to the corresponding value of e. 
One example for such solution is 6 = 1; a = 1/2, e = oo. Then the solution 
of the Lorenz system can be constructed on the basis of elliptic functions. 
We let the interested reader to find the values of parameters for which the 
generalized Lorenz system is integrable. 



Qa^ -ab-2a -b{b-l) = 0, 
-15(6-l)]-9a(2a-l) = 0, 
2 , ..bil-ae^] 



5) I 2aa^ + 6aiCi 
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5 Method of simplest equation and its ver- 
sion called modified method of simplest equa- 
tion 

From the many approaches for obtaining exact analytical solutions of nonlin- 
ear PDEs here we shall give short remarks on the method of simplest equation 
[T6t [T7] and on his version called modified method of simplest equation |18j - 
[20] . The method of simplest equation uses the first step of the test for the 
Painleve property: the leading order analysis to determine the truncation of 
the sum from Eq. (1261) below. In the modified method of simplest equation 
this first step is changed by the equivalent procedure of solving of balance 
equation. The methodology is based on the fact that after application of 
appropriate ansatz large class of NPDEs can be reduced to ODEs of the kind 

( , , dF <PF \ 

(25) p^f«)__..^J=0. 

and for some equations of the kind (125|) particular solutions can be obtained 
which are finite series 

(26) F{0 = jlam)t 

i=l 

constructed by solution $(^) of more simple equation referred to as simplest 
equation. The simplest equation can be the equation of Bernoulli, equation 
of Riccatti, etc. The substitution of Eq. ( 12^ in Eq. leads to the 

polynomial equation 

(27) p = ao + ai$ + a^^^ + ... + a,<l>" = 0, 

where the coefficients cij, i = 0,1,. ..,r depend on the parameters of the 
equation and on the parameters of the solutions. Equating all these coeffi- 
cients to 0, i.e., by setting 

(28) cri = 0,^ = l,2,...,r, 

one obtains a system of nonlinear algebraic equations. Each solution of this 
system leads to a solution of kind (126|) of Eq. (125|) . 

In order to ensure non-trivial solution by the above method we have 
to ensure that cr^ contains at least two terms. To do this in the modified 
simplest equation method we have to balance the highest powers of that 
are obtained from the different terms of the solved equation of kind fl25|) . 
As a result of this we obtain an additional equation between some of the 
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parameters of the equation and the solution. This equation is called balance 
equation [IB]-[2n]. 

As conclusion we stress that there is a unproved conjecture that the 
Painleve property is a direct test of integrability of a differential equation 
or for a system of differential equations. There is however a solid evidence 
that the conjecture is true [12j. The Liouville theorem from the complex 
analysis states that if f{z) is entire bounded function of the complex vari- 
able z the it can take only one (constant) value f{z) = c. Example of such 
entire bounded functions (if we think the time t to be complex variable) 
are the time-independent integrals of the motion connected to a differential 
equation or to a system of such equations. Let us assume that the studied 
system of equations is Hamiltonian and algebraic integrable, i.e., the time- 
independent integrals are polynomials of the canonical variables Pi and g^. In 
the complex plane the functions pi and qi can have various kinds of movable 
singularities, but...., as the constructed by them integral has to be a constant 
in the complex plane all the singularities must cancel at any singularity po- 
sition to- But this is possible for polynomial in pi and qi integrals only if the 
singularities of Pi and q^ are simple poles or rational branch points. For other 
kinds of terms in the integrals it will be not possible to construct a polyno- 
mial function in which the corresponding power of singularities cancel. The 
words above hint that if the studied system has Hamiltonian with integrals 
that are not polynomials (and are for an example irrational or transcendental 
functions) the Painleve property could not be considered as appropriate test 
for integrability of corresponding system of differential equations. 
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